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Abstract 

The mesosocpic concept is applied to the theory of mixtures. The 
aim is to investigate the diffusion phenomenon from a mesoscopic point 
of view. The domain of the field quantities is extended by the set of 
mesoscopic variables, here the velocities of the components. Balance 
equations on this enlarged space are the equations of motion for the 
mesoscopic fields. Moreover, local distribution functions of the veloci- 
ties are introduced as a statistical element, and an equation of motion 
for this distribution function is derived. From this equation of motion 
differential equations for the diffusion fluxes, and also for higher order 
fluxes are obtained. These equations are of balance type, as it is postu- 
lated in Extended Thermodynamics. The resulting evolution equation 
for the diffusion flux generalizes the Fick's law. 

1 Introduction 

Most phenomena involving mass diffusion are described by means of Fick's 
law which relates linearly the mass flow to the gradient of the mass concen- 
tration through the diffusion coefficient. In practical applications Fick's law 
is satisfactory to a large extent. However, it is not applicable to transitory 
situations involving for instance high frequencies. The mass balance together 
with Ficks law leads to parabolic differential equations for the concentra- 
tions, allowing for infinite velocities of propagation of disturbances. The last 
paradox must be overcome by considering non-Fickian transport behavior 
for the diffusive mass flux. A possible solution is given in Extended Ther- 
modynamics 121 El 01 Ej • The diffusion flux is introduced as an additional 
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field variable, obeying a differential equation. This differential equation is 
derived either from an Irreversible Thermodynamics treatment E] of the 
dissipation inequality (for other applications see [HI H3 E3 EEU El ECU ) , or a 
balance type equation for the non-convective fluxes, here the diffusion flux, is 
postulated in Rational Extended Thermodynamics ^3J 03 El ■ For a purely 
macroscopic derivation of balance type differential equations for higher or- 
der fluxes see 117) . Within Rational Extended Thermodynamics it has 
been shown that the resulting system of field equations is symmetric hyper- 
bolic with a convex extension \1H\ 15} ITH]. thereby allowing for finite speeds 
of disturbances only fZ ty TH ] ITS]. 

A foundation of the field equations of Extended Thermodynamics has 
been given from statistical mechanics. In most cases the single particle dis- 
tribution function is used, together with Boltzmann-equation 0. In this 
case microscopic interactions are introduced via the collision term in the 
Boltzmann-equation. However, Boltzmann-equation is, strictly speaking, 
valid only for dilute gases, and for liquids or solids a different foundation 
is necessary. Such a background theory can be given using the BBGKY- 
hierarchy ( |221 12.S1 1241 125j ) for the hierarchy of N-particle distribution func- 
tions. The structure of the field equations of Extended Thermodynamics 
could be derived from this statistical mechanics background. Such an ap- 
proach necessitates the knowledge of microscopic inter-particle interaction 
potentials. 

An approach, different from microscopic statistical mechanics, is the so 
called mesoscopic theory j^Hl HH1 ED] • It is a continuum theory, introduc- 
ing constitutive functions on a continuum level, and no microscopic inter- 
particle interactions. On the other hand the mesoscopic description is finer 
than the macroscopic description. This additional information leads to the 
definition of internal variables in complex media, and to equations of mo- 
tion for them. This mesoscopic concept has been applied to liquid crystals 
[SH E21 EH1 EH ESI EHl EZI , to damage and fracture mechanics [H2EB], and to 
dipolar media [2E] • We will show in the following that within the mesoscopic 
theory a balance type equation for the diffusion flux can be derived. 
The basic idea of the mesoscopic theory is to introduce an additional vari- 
able, the so-called mesoscopic variable, in the domain of the field quantities. 
On this higher dimensional space a continuum theory is developed. In this 
theory the mesoscopic field quantities depend not only on position and time, 
but also on the mesoscopic variable. In addition a local distribution function 
of the mesoscopic variable is introduced as a statistical element. Macroscopic 
fields are the averages of mesoscopic ones over the additional (mesoscopic) 
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variable. 

The aim of this paper is the application of the mesoscopic concept to 
an N-component mixture and a derivation of differential equations for the 
diffusion fluxes from the mesoscopic theory. The application of the resulting 
equations to fast phenomena, involving diffusion processes is left for future 
work, as well as the investigation whether our generalized diffusion equation 
leads to hyperbolic field equations. 



2 Mesoscopic concept and balance equations 

The mesoscopic concept introduces the mesoscopic space 

(m, x,f)eMxR^xKi 

on which all field quantities are defined. Here m is a set of mesoscopic vari- 
ables which is an element of a suitable manifold M. on which an integration 
can be defined. Beyond the use of additional variables m the mesoscopic con- 
cept introduces a statistical element, the so-called mesoscopic distribution 
function (MDF) /(m,x, t) generated by the different values of the meso- 
scopic variable of the particles in a volume element 

/(m, x, t) = /(•), (•) = (m, x,t)eMx Rl x R t . (2.1) 



The MDF describes the distribution of m in a volume element around x at 
time t, and therefore it is normalized 



/ 



/(m,x,t)dm = 1. (2.2) 



It is the probability density of having the particular value m of the meso- 
scopic variable in the continuum element at position x and time t. Therefore 
it is the fraction 

p(m,x,f) 

/(m,x,i) = — ■ 2.3) 

Here p(x, t) is the macroscopic mass density. By use of (|2.2|) we obtain 

p(x,i) = J p(m,x,i)dm. (2.4) 



4 



2.1 Balance equations 

Let G denote a region in M 3 x A4 and X the density of an extensive quantity. 
Then the global quantity in the region G changes due to a flux over the 
boundary of G and due to production and supply within G: 

— [ Xd 3 xdm= [ (f> x (-)da+ [ E x (-)d 3 xdm . (2.5) 
dt Jg JdG Jg 

A generalized Reynolds transport theorem in the mesoscopic space, analo- 
gous to the one in is used to transform the time derivative, and Gauss 
theorem is applied to the boundary integral. The boundary dG of G consists 
of a boundary in position space and a boundary on the manifold A4. Then 
we have in regular points of the continuum the local mesoscopic balance [IU] 

|x(.) + V, • [v(.)X(.) - 5(-)] + V m • M-)X(-) - R(-)] = £(•)• (2-6) 

Here the independent field w (•), defined on the mesoscopic space, describes 
the change in time of the set of mesoscopic variables: With respect to m the 
mesoscopic change velocity w(-) is the analogue to the mesoscopic material 
velocity v(-) referring to x: If a particle is characterized by (m,x,t), then 
for At — > +0 it is characterized by (m + w(-) At, x + v(-) At, t + At) at later 
times t + At. Besides the usual gradient, also the gradient with respect to 
the set of mesoscopic variables appears. This expresses the fact that there 
is a flux over the boundary on the manifold Ai. 

According to equation (J2.2|) we obtain from the mesoscopic mass balance a 
balance type differential equation for the MDF /(•) by inserting its definition 

duinni: 

^/(•)+V»-[v(.)/(-)]+V m >(•)/(•)] + 



+/(• 



a + v(-)-V. 



]ng(x,t) = 0. (2.7) 



3 Application of the mesoscopic concept to mix- 
tures 

Let us consider a mixture of N different chemical components. The different 
components are distinguished by a component index in capital letters. We 
will introduce velocity distributions of all chemical components, denoted f A . 
Each of them is a function of position of the respective continuum element, 
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time, and particle velocity. f A (x, t, v) gives the probability density of finding 
the value v for the material velocity of component A of the mixture in a 
volume element around x at time t. 

The mesoscopic mass densities and all other mesoscopic fields are differ- 
ent for different chemical components. Mesoscopic fields are denoted by a 
" hat" , where this is necessary in order to distinguish them from the corre- 
sponding macroscopic quantities. 
The distribution function f A is defined as follows: 



/ (x, t, v)dv gives the fraction of particles of component A having a material 
velocity in the element dw on ft^ around v. 
The distribution functions are normalized: 



4 Macroscopic balance equations for the mixture 
and its components 

In this section, we shall recall briefly the macroscopic balance equations for 
mixtures and for the different components. 

The total mass of a mixture is a conserved quantity, while the mass of 
a particular component is not conserved, if chemical reactions occur. The 
same happens for momentum, and energy due to chemical reactions between 
different components of the mixture. Therefore, in the balance equations for 
components there are production terms on the right hand side. 

4.1 Macroscopic balance equations for the mixture 

The balance equations for the mixture as a whole look the same as the 
balance equations for a one-component system |421 I43| . Prom the struc- 
ture of balance laws, mixtures cannot be distinguished from chemically pure 
substances. 

Balance of mass 




(3.1) 




(3.2) 




(4.1) 
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Here v is the macroscopic material velocity of the mixture, i.e., the barycen- 
tric velocity of all particles in the volume element. 
Balance of momentum 

— (pv) + V-( / 9v®v-t T ) = pf , (4.2) 

where t denotes Cauchy stress tensor, and f is the acceleration (specific den- 
sity of volume forces). Balance of energy 



d 1 1 

+ -v • v)] + V • [q - v • t + pv(e + -v • v)] 

-p(f-v + r)=0 (4.3) 

where e is the specific internal energy density (internal energy per unit mass), 
q denotes the heat flux density, and r is the energy supply density. 



4.2 Macroscopic balance equations for the components 
Balance of mass of component A 

The local balance of mass of component A, taking into account production 
of mass of that component due to chemical reactions, reads: 

^(x,t) + V ■ (/v^)(x,t) = P c A hem (4.4) 

here v A is the material velocity of component A. The mass densities of the 
different chemical components are additive: 



£/ = P (4-5) 



The material velocity v of the mixture is defined as the weighted sum of the 
component velocities: 

Y,P A v A = pv (4.6) 

A 

Equation (|4.4j) may be cast in a more useful form in terms of mass fractions 
c A defined by c A = A = 1, N and the diffusion fluxes 3 A 

J A = p A (v A -v) (4.7) 

such that Y.A jA = °- The 

mass fraction balance equation is 

Pj t c A + V-J A = P A hem (4.8) 
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where ^ = Jj + v • V is the material time derivative. 
Balance of momentum of component A 

The momentum balance for one particular chemical component is: 

^(p A v A ) + V • (p A v A ® v A - t AT ) = pH A + V A . (4.9) 

with a production P A of momentum of component A. 
Balance of energy of component A 

The energy balance equation is given by: 

^[p A (e A + ^v A • v A )] + V • [q A - v A • t A + p A v A (e A + K A • v A )] 
-/(f^-v A + r A ) =P A (4.10) 

with a component energy production due to chemical reactions. In this 
set of balance equations constitutive equations for the component stress 
tensor, component internal energy density, and heat flux density, and the 
productions due to chemical reactions are needed in order to close the system 
of differential equations. These constitutive quantities must be related to 
the wanted fields in a material dependent manner. The constitutive theory, 
including restrictions on constitutive functions from general principles, is 
out of the scope of the present paper. 

The (total) energy density of the mixture is the sum over the compo- 
nents. From the additivity of the extensive quantities mass, momentum, and 
energy follows that the balances of the mixture are obtained summing up the 
component equations [JSE!]- This leads to relations between constitutive 
quantities for components and constitutive quantities for the mixture. 

5 Mesoscopic balance equations for the mixture 
and its components 

Mesoscopic balance equations are defined on the mesoscopic space. For the 
mixture as a whole the additional mesoscopic variable is the material velocity 
v of the mixture, which is the weighted sum of the component velocities, see 
equation Q4.6|) . The mesoscopic balance equations are derived analogously 
to the derivation discussed in the second section (see equation (|2.6j) for the 
general structure of these balances) and by suitable identifications. The 
mesoscopic balance equations for the mixture are those derived previously 
for one component systems (see Section 2). In all cases it is supposed that 
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the constituents of the mixture have no internal angular momentum, i.e., 
we are not dealing with micropolar continua. 

For the mixture we introduce the distribution function of the material 
velocity: 

Explicitly we have the following mesoscopic balances: 



5.1 Mesoscopic balance equations for the mixture 
Mesoscopic balance of mass 

^(P)(«) + V • (pv)(.) + V v • (wp)(.) = (5.2) 

here (•) is the abbreviation for (x,t,v), v and w are the mesoscopic ve- 
locity of the mixture and the change in time of the mesoscopic variable v, 
respectively (w = v). 

According to the definition of the distribution function equation (|5.1|) . we 
obtain from the mesoscopic mass balance (|5.2j) a differential equation for 
the MDF /(•): 

?M + V • (/(.) V ) + V v (/(.)w(.)) + (5.3) 
+/(-)(^ + v-v) lgp(x,i) = 
Mesoscopic balance of momentum 

— (pv) + V • (pv <g) v) + V v • (pw ® v) - V • f - V v • T T = pf. (5.4) 

with t = t and T = T. Here T is the analogue in velocity space of the 
stress tensor, i.e., the non-convective momentum flux in velocity space. All 
field quantities here depend on position, time and particle velocity. 

Due to definition ()5.1|) the macroscopic material velocity of the mixture, 
v(x, t), i.e., the mean velocity of the mixture, is expressed as follows: 

p(x,t)v(x,f) = / pvdv, (5.5) 
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Mesoscopic balance of energy 

d 1 .1 

g- t [p{e + -v • v)] + V • [q - v t + p(e + -v • v)v] + 

+V V ■ [Q - w T + pw(e + • v)] = 

= p(f ■ v + f) + P e (5.6) 

Q is the analogue to the heat flux density in velocity space, i.e., a non- 
convective energy flux over the boundary in the mesoscopic part of the 
domain, which consists of position, time and velocity for all fields, here. 

By integrating these mesoscopic balances over the mesoscopic space, i.e., 
over all velocities, with the meosocopic distribution function as statistical 
weight, the macroscopic balance equations for the mixture are obtained. 



5.2 Mesoscopic balance equations for the components of the 
mixture 

Mesoscopic balance of mass of component A 

By using the general form of the mesoscopic balance of mass (see also |331 
1441 1411 IHUj ^. we have the following equation for the component A: 

|(p A )(.) + V • (/v A )(-) + V v • (wV)f') = PLrai-) (5-7) 

here v A represents the material velocity of component A, and w A is the 
change in time of the mesoscopic variable v A : w A = v A . This acceleration 
w A is a mesoscopic constitutive function. The term on the right hand side 
represents a production due to chemical reactions. 

The sum of the mass balances of the different components is the mass 
balance of the mixture, because the mass densities are additive: 

E p A = p- ( 5 - 8 ) 

A=1,...,N 

I E(^) + v • B A A ) + v v • J> V) = E p ctm ■ (5-9) 

A A A A 

By comparing in equations Q5.9|) and (|5.2j) the fluxes in mesoscopic space 
and in position space we obtain the following relations: 

E V v • (w V) = V v • (wp) , (5.10) 

A 

^ V .(p A v A ) = V-(pv) , (5.11) 

A 
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and for the production terms: 




' j *-chem 
A 



(5.12) 



A solution of equations Q5.1U|) and (|5.11|) is: 



w = 



v 




(5.13) 



(5.14) 



Mesoscopic distribution function for the components 

The distribution function f A (x, t, v) gives the probability density of finding 
a particle of component A with velocity v in the volume element around 
x at time t. It is the mass fraction 



The mesoscopic balance of mass for component A, equation (|5.7jl . can be 
used to derive an equation of motion for the distribution function: 



Introducing the macroscopic fields by integrating over the mesoscopic part 
of the enlarged domain, here the velocity, with the component distribution 
function as probability density, we can transform local mesoscopic balances 
into local macroscopic balances. So, the definition of the macroscopic ma- 
terial velocity of the component A, v^(x, t), i.e., the mean velocity of the 
component A of the mixture, is 



and equation ()5.7j) leads to the macroscopic mass balance of component A 
if we take 




(5.15) 




(5.16) 




(5.17) 




(5.18) 
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into consideration. This holds because of Gauss' theorem and because / (•)> 
restricted to the mesoscopic part has a compact support. 



Relation between the mesoscopic distribution functions of the 
components and the distribution function of the mixture 

The velocity distribution function of a particular component / was 
defined as 

and the velocity distribution of the mixture is 

/W,v, = ^. (5.20) 

Because the mass densities of the different chemical components are ad- 
ditive: 

/(x, t, v) = p(x, t, v) , *) = ^( x ' *) > ( 5 - 21 ) 



we end up with a relation between the distribution functions of the compo- 
nents and that of the mixture: 

/ A (x, t, v)/(x, t)=J2 p A (x, t, v) = p(x, i, v) = /(x, i, v)p(x 5 1) , 
A A 

(5.22) 

or 

£/ A (x,t,v)^^=/(x,t,v). (5.23) 

The distribution function of the mixture is the sum of the component dis- 
tribution functions, weighted with the mass fractions of the components. 
Mesoscopic balance of momentum of component A 

^(^V)+V'(A A ®v A )+Vv(A A ®v A )-V.t A -V v -f A = p A F A +t>i. 

(5.24) 

The mesoscopic momenta p A v A are additive, and summing up equations 
(|5.24|) over the different components one obtains the mesoscopic balance of 
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momentum (|5.4|) . By comparing the fluxes in position space and in velocity 
space, respectively, we have: 

t = ^(i A - p A 5v A ® <5v A ), (5.25) 

A 

T = ^(T A - / S A 5w A ®5v A ), (5.26) 

A 

where we have introduced the abbreviations 

5v A = v A - v, 5w A = w A - w. (5.27) 

Finally, a comparison of the production terms in the mesoscopic balances of 
momentum for the sum of components on one hand, and for the mixture on 
the other hand, leads to the relation: 

pf = J2(P A f A + K) (5-28) 
A 

By integrating over the mesoscopic space the macroscopic momentum bal- 
ance for component A, equation ()4.9|) is obtained. 
Mesoscopic balance of energy of component A 
For particles of component A and particle velocity v A we have: 

|[/(e A + \v A ■ v A )] + V • [q A - v A • i A + p A (e A + \v A • v> A ] + 

+V v a • [Q A - w A • f A + p A ~w A (e A + \v A ■ v A )] = 

= p A (i A ■ v A + f A ) + P A (5.29) 

As before, summing over all components of the mixture one obtains the 
mesoscopic balance of energy (|5.6|) and the integration over IR 3 leads to the 
macroscopic balance of energy expressed by equation (|4.3|) . 

6 Constitutive quantities for the components com- 
pared to those for the mixture 

Analogously to the previous section the flux terms and the production terms 
in the balance equations of the mesoscopic mixture and the sum of the 
component equations are compared. As a result the internal energy density 
of the mixture, the heat flux and its analogue in mesoscopic space, and 
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the energy absorption density of the mixture are expressed by the following 
relations: 

pe = Y,{p A Z A + \p A 8v A .8v% (6.1) 

A Z 

q = ]T{q A - 5w A ■ t A + p A e A 8w A + \[p A {5v A ■ 5v A )5v A ]}, (6.2) 
A 2 

Q = ^2{Q A - p A (w • 5w A )5v A - 5w A ■ T A + p A e A 6w A + 

A 

+p A (5v A • v)<5w A + ~[p A (8v A ■ 5v A )5w A ]}, (6.3) 

pr = J2(f> A r A + p A v A ■ 8f A ) with 8f A = f A - f , (6.4) 
A 

A 

Equations (|6.1|) to Q6.5|) show that the mesoscopic constitutive quantities of 
the mixture are not, in general, the sum of the corresponding constitutive 
quantities of the components, but some fluctuation terms contribute. 



7 Evolution equations for diffusion fluxes 

The aim of this section is to obtain evolution equations for diffusion fluxes of 
the different components. By using relations (|5.5|) . Q5.17|) and the definitions 
(|3.1|) . ()5.1I) it is possible to define the diffusion flux J A , expressed by equation 
(|4~7I) . as follows: 

J A (x,t) = /(x,t) / [f A (-)-f(-)]vdM (7.1) 

In order to derive evolution equations for fluxes let us introduce the family 
of the macroscopic fields of order parameters which is defined by different 
moments of the distribution functions f A (-) and /(•). They are defined as 

4=! f A (-)V7^dM, (7.2) 

a k = f f(-)V7TVdM, (7.3) 
Jm.3 — * — 
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where nrr 1 denotes the symmetric irreducible (traceless) part of a tensor 
|33| . These fields of order parameters describe macroscopically the meso- 
scopic state of the system introduced by v A and v and its distribution 
functions f A and /. Thus these are the link between the mesoscopic back- 
ground description of the system and its extended description by additional 
macroscopic fields. In a macroscopic phenomenological theory they repre- 
sent in general internal variables which satisfy relaxation equations. In our 
case the first moments ai and a A are not internal variables, but classical 
wanted fields, namely the specific momentum densities (material velocities) 
v and v A . The higher order moments are internal variables in the sense of 
thermodynamics. 

The diffusion fluxes 3 A , (A = 1, N) are proportional to the difference 
of the first moment of the distribution functions / (■) and /(■). If we mul- 
tiply equations (|5.16j) and ()5.3|) with the mesoscopic variable and integrate 
over the manifold A4 we obtain, respectively: 



/ ^[f A (-)v A ]dM+ [ V-[v A f A (.)v A }dM + 

Jm.3 ot J R3 



+ / v A V v • [w A (-)f A (-)]dM + / / A v 

J UP 

pA 
chem _ dM 



dt 



+ v • V 



lg p A (x,t)dM 



p A (x,t) 



(7.4) 



and 



i/(.)v]dM+ / V-[v/(-)v]dM 
a at j R 3 



+ / vV v • [w(-)f(-)}dM + / /v 
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^^(x,*)^ = 0. 



The mesoscopic manifold M, here the M 3 spanned by the velocities v, is 
time independent. The time derivative and the derivative with respect to 
position can be interchanged with the integration over A4. We split the 
mesoscopic velocity into the macroscopic velocity and the deviation from 
this average, respectively, for component A and the mixture: 



v A (x,t) +5v y 



(7.5) 



and 

v(-) = v(x,t) + «5v. 



(7.6) 



15 



The resulting equations of motion for the first moments a A and ai are, 
respectively: 



^af + V- [v A (x,i)a^+ / 5v A (-)f A (-)v A dM 
at - j R 3 

d 



+ 



(7.7) 



v^V v • [f A (-)w A (-)]dM + af [- + v A (x, t) ■ Vj lg /(x, t) + 



Vlg/(x,t) / 5v A (-)f A (-)v A dM 



1 



P A (x,t) 



P C A hemV A dM 



and 



v(x,t)ai+ / (5v(-)/(-)vdM 



3 



vV v • [/(-)w(-)]dM + a x [- + v(x, t) ■ Vj lg p(x, t) + (7.f 



+Vlgp(x,t) / <yv(.)/(-)vdM = 



We now introduce some approximations: the deviations 6v A and <5v vanish, 
i.e., the mesoscopic variable v has the value of the barycentric velocity. In 
other words we suppose v (•) = v j4 (x, t) and v(-) = v(x, t). 
By computing the difference of equations <|7.7jl and (|7.8|) and taking into 
account these approximations one has: 



d f 
— (& A - ai) + V • (v A a A - vai) + / vV, 
ct Jm 3 

-/(.)w(.)dM + a^[|- + v A (x, i) • vl lg /(x, *) 



ai 



.at 

| + v(x,t).v]lg,(x,t) = -^ 



P, 



chem 



dM 



(7.9) 



The last equation (j7.9j) represents a general phenomenological evolution 
equation for diffusion fluxes in an "extended theory". Using further ap- 
proximations it is possible to obtain a Cattaneo-type diffusion equation. In 
order to exploit equation (|7.9|) further, it is necessary to insert expressions 
for w and w A . These are equations on the mesoscopic level, and they can 
be interpreted as mesoscopic constitutive equations. In this case the do- 
main of constitutive mappings, the state space has to be introduced. In 
mesoscopic theories there are different possibilities: one is to introduce a 
state space of purely mesoscopic variables, another one is to introduce a 
mixed space including mesoscopic and macroscopic variables, and the third 
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possibility is one containing only macroscopic variables. The mesoscopic 
background of the constitutive theory has been applied to liquid crystals of 
uniaxial molecules j3U S3 IS] and to dipolar media [IT]- For diffusion 
in mixtures this is left for future work. 

8 Conclusions 

We have shown a way to derive a differential equation for the diffusion flux 
from the so called mesoscopic theory. In this refined theory an additional 
variable in the domain of the field quantities is introduced, here the velocity, 
and we have a velocity distribution. We have derived a differential equation 
for the diffusion flux from the mesocopic balance equations. The differential 
equation for the diffusion flux is of balance type, i.e., of the form used in 
extended thermodynamics. Therefore, it can be expected that our approach 
leads to hyperbolic systems of field equations, but the investigation of this 
question is left for future work. 

Because the level here is not the microscopic one, no assumptions about 
inter-particle interactions have to be made here in contrast to kinetic theory. 
We are here on the continuum level, and many-particle interactions are taken 
into account automatically. This is not the case for derivations of the field 
equations of extended thermodynamics based on the Boltzmann-equation. 
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